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Abstract. 

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) 
and its interaction with the coronal magnetic field based on the May 12, 1997, CME event 
using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambi- 
ent solar wind steady-state solution is driven by photospheric magnetic field data, while 
the solar eruption is obtained by superimposing an unstable flux rope onto the steady- 
state solution. During the initial stage of CME expansion, the core flux rope reconnects 
with the neighboring field, which facilitates lateral expansion of the CME footprint in 
the low corona. The flux rope fleld also reconnects with the oppositely orientated over- 
lying magnetic field in the manner of the breakout model. During this stage of the erup- 
tion, the simulated CME rotates counter-clockwise to achieve an orientation that is in 
agreement with the interplanetary flux rope observed at 1 AU. A significant component 
of the CME that expands into interplanetary space comprises one of the side lobes cre- 
ated mainly as a result of reconnection with the overlying field. Within 3 hours, recon- 
nection effectively modifies the CME connectivity from the initial condition where both 
footpoints are rooted in the active region to a situation where one footpoint is displaced 
into the quiet Sun, at a significant distance (w ^Rq) from the original source region. 
The expansion and rotation due to interaction with the overlying magnetic field stops 
when the CME reaches the outer edge of the helmet streamer belt, where the field is 
organized on a global scale. The simulation thus offers a new view of the role reconnec- 
tion plays in rotating a CME flux rope and transporting its footpoints while preserv- 
ing its core structure. 
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1. INTRODUCTION 

In recent decades, the concept of magnetic reconnection, 
in which two oppositely orientated magnetic field hnes are 
cut and re-assembled, has become more accepted as an im- 
portant process in solar coronal dynamics. In particular, 
magnetic reconnection seems to dominate the initiation, evo- 
lution, and propagation of Coronal Mass Ejections (CMEs), 
as the magnetic flux carried by the CME interacts with the 
pre-existing field both in the low corona and throughout the 
heliosphere. In recent years, a number of CME initiation 
models have attempted to resolve the observed features of 
a solar eruption. Many CME initiation models drive the 
eruption by imposing photospheric shear motion, flux dif- 
fusion, and flux cancellation to the footpoints of the pre- 
existing magnetic field (see [e.g. Forbes and Priest, 1995; 
Amari et al, 2000; Linker et ai, 2003] and the review by 
Forbes et al. [2006] and references therein). In other words, 
the eruption is obtained through modiflcation of the bound- 
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ary conditions, transforming the initial potential magnetic 
field into a non-potential, twisted magnetic field, which has 
the required free energy stored in it. This non-potential 
state introduces currents so that the Lorentz force over- 
comes the downward magnetic tension and gravity; eventu- 
ally, the CME erupts when the force balance breaks down. 
In this class of CME initiation model, which we refer to 
as the "driven" model, magnetic reconnection occurs below 
the CME flux rope in the vicinity of the magnetic field foot- 
points. 

Another type of CME initiation model is the breakout 
model [Antiochos et at, 1999; Lynch et al, 2008], in which 
a bipolar region (representing the source active region of the 
CME) is embedded in a background dipole fleld of opposite 
polarity. The eruption in this case is driven by shearing 
the footpoints of the local bipolar region around its neu- 
tral line. As a result, the local bipolar field inflates and 
begins to reconnect with the oppositely orientated overlying 
field. This reconnection transfers overlying fleld to the side 
lobes, until the CME erupts through the weakened overlying 
field. The location of magnetic reconnection that drives this 
type of CME initiation model occurs above the CME flux 
rope, and there is no opening or disconnection of flux dur- 
ing the process. There are two main problems with most of 
the current CME initiation models. First, these models use 
idealized magnetic conflgurations, which only mimic the re- 
alistic coronal fleld, along with idealized footpoint motions. 
Second, almost all models neglect the non-potential, back- 
ground solar wind solution into which the CME is erupting 
[Manchester et al, 2004; Lugaz et al., 2007; Cohen et al., 
2008a]. Taking into account a more realistic coronal environ- 
ment might significantly affect the propagation of the CME 
due to the interaction with the more complex ambient field. 
This three-dimensional interaction can facilitate the lateral 
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expansion of the outer shell of the CME via magnetic re- 
connection during the eruption [e.g. Manoharan et al., 1996; 
Pick et al, 1998; Pohjolamen et al, 2001]. 

One of the observed signatures associated with CMEs 
in the low corona are the so-called coronal "waves". Since 
their discovery in 1996 [Dere et al, 1997: Moses et al, 1997; 
Thompson et al, 1998], the physical nature of EIT coro- 
nal waves, which are strongly linked to CMEs [Biesecker 
et al, 2002] has been under debate. In particular, there is 
an argument whether the observed coronal "wave" is a Mag- 
netoHydroDynamic (MHD) wave, or whether it corresponds 
to the actual footprint of the CME (for a complete discus- 
sion on this debate see Cohen et al [2009], as well as Ofman 
and Thompson [2002]; Ofman [2007]; Sdvmidt and Ofman 
[2010] and references therein). In one non-wave model, the 
latter requires the flux rope to reconnect with the surround- 
ing coronal field, facilitating the lateral expansion [Attrill 
et al, 2007]. This paper demonstrates that the concept can 
occur in conjunction with the breakout model where the 
original flux rope reconnects with the overlying fleld. The 
main difference between the two is essentially the location of 
the recormcction point. In the breakout model, the polarity 
of the CME flux rope is opposite to that of the large-scale 
overlying field, so that the reconnection occurs at the top of 
the flux rope. In the model by Attrill et al [2007], which 
we can call the "stepping reconnection model" , reconnection 
occurs whenever the core flux rope meets a neighboring loop 
with opposite polarity and the reconnection point is located 
to the side of the expanding core flux rope. We demonstrate 
that the stepping reconnection model can theoretically oc- 
cur within the central part of the breakout model topology, 
facilitating the expansion of the core flux rope, prior to in- 
teraction with the larger-scale overlying field. This concept 
is discussed in Section 4. 

In a three-dimensional MHD simulation of a recent CME 
event (February 13, 2009), Cohen et al. [2009] have shown 
that the core flux rope indeed reconnects with the surround- 
ing field in the same manner as described by Attrill et al. 
[2007]. They also showed that the bright front constituting 
the diffuse EIT coronal waves is composed of both a piston- 
driven MHD wave as well as non-wave components, which 
are coupled as long as the CME continues to expand later- 
ally. The lateral expansion essentially stops when the field 
topology no longer allows the magnetic field to reconnect, 
and the magnetic overpressure of the flux rope relative to 
the surroundings no longer drives a strong lateral expansion. 
From this point, which can be at the considerable distance 
of w IRq from the source region [Attrill et al, 2009; Cohen 
et al, 2009], only the MHD wave component continues to 
exist as a freely propagating wave. 

The May 1997 CME event is a SHINE campaign event 
and it has been chosen due to the fact that it was an isolated 
halo CME event occurred during solar minimum, which sup- 
posedly makes it a "simple" event, and due to the relatively 
extensive observational data of the event. Several attempts 
to simulate the global evolution of a CME event in the 
corona and in the heliosphere have been done in the past 
decade. For example, Manchester et al. [2004]; Fan and Gib- 
son [2007]; Rtley et al [2008]; van der Hoist et al [2009] (and 
references therein) have sirrmlatcd CME eruptions driven by 
different methods into cither a background corona in hydro- 
static equilibrium or a more physical solar wind background. 
However, even the steady state solar wind solution in these 
simulations was based on an idealized dipole configuration 
for the solar magnetic field. The particular May 1997 event 
has been simulated by Odstrcil et al [2004, 2005]; Shen et al. 
[2007]; Wu et al [2007]. All these simulations however, were 
focused on the propagation of the CME to the Earth, and 
they were driven by kinematic propagation of the MHD pa- 
rameters to the inner boundary of the simulation domain, 
which has been set to be beyond the Alfvenic point. 

In this work, we use a model for the solar corona which is 
driven by high-resolution magnetogram data and provides a 



more realistic coronal magnetic field. The model also sim- 
ulates a CME that propagates through a steady-state, non- 
potential MHD solution for the solax corona and the solar 
wind. The advantage of such a model is that it enables us to 
study the complex interaction of the CME with the realistic 
ambient field. The main limitation of the model is that the 
CME is a fully-formed fiux rope that is "injected" into the 
steady-state solution with its observed parameters. While 
reconnection between the closed loops of the fiux rope and 
the surrounding magnetic field is a primary focus of this 
paper, any reconnection between closed loops beneath the 
CME that generates flux rope coils, as occurs in many CME 
initiation models, is not part of our simulation. Another dif- 
ference between this model and the CME initiation models 
described previously is that here we are interested in the in- 
teraction above the surface with fixed photospheric bound- 
ary conditions, while in the other models, modification of 
the boundary conditions is the main driver for the sinmla- 
tion. We present a high-resolution numerical simulation of 
the complex interaction between the erupting CME and the 
ambient coronal field (up to 24J?©) based on the May 12, 
1997, CME event. The goal of this work is to better under- 
stand the three-dimensional interaction of the CME with 
the ambient flux. 

The structure of this paper is as follows. We describe the 
numerical simulation in Section 2 and present the results in 
Section 3. We discuss the consequences and implications of 
the simulation results in Section 4 and conclude our findings 
in Section 5. 

2. Numerical IVIodel 

2.1. Ambient Solar Wind and CME Initiation 

For our study, we use the Solar Corona (SC) module [Co- 
hen et al, 2007] of the Space Weather Modeling Framework 
(SWMF) [Toth et al, 2005], which is based on the BATS- 
R-US global MHD code [Powell et al, 1999]. The steady- 
state solar wind solution is obtained in a non-polytropic 
manner [Roussev et al, 2003b; Cohen et al, 2007, 2008b] 
and it is constrained by the empirical Wang-Sheeley-Arge 
(WSA) model [VFanp and Sheeley, 1990; Arge and Pizzo, 
2000). We use a potential magnetic field [Altschuler and 
Newkirk, 1969] to prescribe the initial condition for the 
magnetic field using high-resolution MDI magnetograms 
(http://sun.stanford.edu/synop/), and the steady-state is 
obtained by iterating the MHD equations: 

If + V • (pu) = 0, 

Pt + v-(puu + p/+£j-sa)=pg, 

^+V-(uB-Bu) = 0, (1) 

V • (ipn^u + ^pu + = Pis ■ u), 

where p, u, B, p, g, and 7 are the mass density, velocity, 
magnetic field, pressure, gravity, and the polytropic index, 
respectively, until convergence is achieved. 

We initiate the CME by superimposing an unstable, semi- 
circular fiux rope based on the analytical model by Titov 
and DernouUn [1999] on top of the ambient solution [Rous- 
sev et al, 2003a]. The fiux rope properties are matched to 
fit the observed properties of the source active region and its 
inversion line. The free energy is controlled by an additional 
toroidal field that produces the observed linear speed of the 
CME. The simulation here is based on a previous simulation 
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done by Cohen et al. [2008a], which includes full propaga- 
tion of the CME to 1 AU and comparison with in-situ data. 
The CME initiation method described here has been used to 
study processes in the solar corona [e.g. Lugaz et al., 2007; 
Manchester et al, 2008]. 

2.2. Simulation Setup 

In order to follow the topology and evolution of the mag- 
netic interaction between the CME and the ambient field, 
we design the grid with high resolution around the active 
region and along the line of CME propagation. The grid 
size around the active region is 1Q~ Rq, and the grid size 
along the propagation line is 1O~^7?0. The grid resolution 
is coarser in other regions of the simulation domain, which 
are not of interest in the context of this work. The time 
step in the simulation is determined by the Alfven speed 
and the grid resolution. The higher the magnetic field and 
smaller the grid size, the smaller the time step. In order 
to compensate for the extremely small time steps in active 
regions, where magnetic fields arc very strong, we set the in- 
ner boundary of the sirrmlation to be at a height of Q.QQRq 
above the surface, where the magnetic field is weaker. 

Figure 1 displays the initial stage of the simulation. The 
left panel shows the large-scale structure of the steady-state 
corona, with selected closed field lines in blue and open field 
lines in red. The background color contours represent the 
solar wind radial speed, ranging from fast in yellow to slow 
in blue. The right panel shows the vicinity of the active re- 
gion at the initial state of the simulation, where its approx- 
imated location on the Sun is marked by the black square. 
Color contours represent the magnitude of the radial field 
on a sphere at a height of 1.06/?©, red streamlines repre- 
sent three-dimensional magnetic field lines of the superim- 
posed flux rope, and solid white lines mark the grid structure 
around the flux rope. The flux rope itself is represented by 
an iso surface of mass density, p = lO"^"* g cm~^, which is 
greater than the surrounding density at the same height. 

2.3. The Concept of Numerical Reconnection in the 
Global Model for the Solar Corona 

The model solves the set of ideal MHD equations, so no 
physical resistivity is introduced in the equations. There- 
fore, the rate of numerical magnetic reconnection in the 
simulation is, in principle, controlled by numerical diffusion 

term, which is designed to stabilize the numerical solution 
and is not directly related to the physical resistivity of the 
system. The value of the numerical diffusion is proportional 
to the size of the grid cell and the time step via the ra- 
tio Ax'^ /At. In the case of the numerical resistivity in the 
induction equation (where it can affect magnetic reconnec- 
tions) this ratio is proportional to ri/no, where jj,o is the per- 
meability of free space, and r] is the electric resistivity. In 
the solar corona, the typical value for the Lundquist num- 
ber, S = (xoLva/v is 10^^ — lO^"* [Boyd and Sanderson, 
2003], where L is the typical length scale, and va is the 
Alfven speed. The typical value of Aa; is a fraction of a so- 
lar radius and the time step can range between 0.1 — 10 s. 
Therefore, using the same typical values for L and va as 
above, we find that 1 < Sn ~ LvAAt/A'^x <IC S, where 
Sn is the Lundquist number calculated using the rmmeri- 
cal resistivity. This means that the model might result in 
ovcr-rcconncction of the magnetic field in regions where two 
opposite magnetic field lines are pushed towards each other. 
An example for such numerical behavior is the generation of 
"U-shape" detached field lines around the heliospheric cur- 
rent sheet. This issue is resolved by implementing the Roe 
solver in the numerical model, which is a more precise, less 
diff'usive numerical scheme, and is equivalent to the addition 
of one more level of grid refinement (regardless of the actual 
smallest grid size). The Roe solver and its implementation 



to the MHD model are discussed in details in Sokolov et al. 
[2008]. 

When discussing the magnetic reconnection of a CME 
with the coronal ambient field, one should keep in mind that 
this reconnection has a global, macroscopic sense, and that 
the time scale for such reconnection (which can be of the 
same order as the Alfvenic time scale) depends on the CME 
size and speed. In the simulation presented here, we try to 
capture this global interaction between the newly injected 
magnetic flux carried by the CME and the pre-existing field. 
This concept is different from the microscopic description of 
magnetic reconnection that obviously needs a better treat- 
ment than a global MHD model. In this simulation, we focus 
our study on the re-distribution of the global fleld as the new 
flux is introduced into the system, and on how the interac- 
tion between the fields affects the propagation of the CME 
through the corona. We compare our results with observa- 
tions of large-scale signatures for such interaction between 
the fields in the corona. We emphasis that any further men- 
tioning in the text for the term "magnetic reconnection" 
refers to a numerical reconnection in the simulation. 

We simulated the first 20 hours of propagation using the 
PLEIADES cluster at NASA's NAS center. The CME front 
left the simulation domain after approximately 14 hours. 

3. Results 

We separate the analysis of the results into two parts. 
First, we validate the timing and the overall structure of 
the CME with some observations, and, second, we analyze 
the results assuming the overall interaction is satisfactorily 
captured by the simulation. Figure 2 compares the ob- 
served and simulated coronal dimmings and the observed 
and simulated coronal wave front associated with the den- 
sity changes due to the CME expansion and propagation 
through the corona [after Cohen et al., 2009]. The four 
rows display comparisons for the early stages of the sim- 
ulation, 10, 30, 40, and 60 minutes after the eruption on- 
set, from top to bottom. The first column shows SOHO 
EIT 195yi base-difference images, constructed with data 
from http:/ /sohowww. nascom.nasa.gov. The second column 
shows the corresponding synthetic EIT 195j4 base-difference 
images produced by the simulation. These images are the 
line of sight integral: J n'iR{ne,T)dl, where Ue is the elec- 
tron number density and R is the response function for the 
EIT 195^ filter based on spectral synthesis from the CHI- 
ANTI code [Landi et al, 2002]. The complete description of 
the synthetic EIT images in SWMF can be found in Downs 
et al. [2010]. The third column shows base-difference im- 
ages of the simulated mass density on a sphere at height 
of r — I.IRq. The differences are normalized to the pre- 
emption mass density distribution in a similar manner as in 
Cohen et al. [2009]. The last column shows a display that is 
similar to the third column but with the addition of selected 
magnetic field lines. Those in the core flux rope are marked 
in red, those with one footpoint rooted in the original source 
region are marked in blue, and overlying field lines with both 
footpoints located far from the source region are marked in 
yellow. The approximate location of the leading edge of the 
coronal wave disturbance is marked by a white circle in the 
first and the third columns. 

Figure 2 shows that the strongest brightening in the syn- 
thetic EIT images (second column) lags behind the matching 
observed (first column) and sirrmlated (third column) coro- 
nal wave fronts, and it remains there in a location compara- 
ble to the location of persistent brightenings observed in the 
base difference images in the first column [see Attrill et al., 
2007; Delannee, 2009]. Although the locations of the ob- 
served and simulated wave fronts match, the latter appears 
only as a weak density enhancement. On the other hand, 
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the source region of the event appears bright both in the 
synthetic BIT images and in the density difference images. 
A discrepancy between the observed and sirrmlated patterns 
is tfie lack of any rotation of the coronal wave in the sirrmla- 
tion. Podladchikova and Berghmans [2005] and Attnll et al. 
[2007] independently showed that in the EIT observations 
the coronal wave undergoes an overall counter-clockwise ro- 
tation. Another discrepancy concerns the coronal dimmings. 
The NW-SE orientation of the main pair of dimmings in the 
simulated images appears as the mirror reflection of those in 
the observations (NE-S W) . Coronal dimmings located on ei- 
ther side of a post-eruption arcade are understood to be the 
footpoints of the flux rope [e.g. Webb et al., 2000). A possi- 
ble reason for these discrepancies is that the flux rope in the 
simulation is idealized and undergoes only a straightforward 
expansion. Real flux ropes can be much more complicated, 
with twist and writhe that might lead to the observed ro- 
tation of the coronal wave and a different representation of 
the dimmings. 

The global field topology as shown in the last column of 
Figure 2 is similar to the topology extrapolated using a po- 
tential field method Delannee [2009] (her Figure 2). This is 
not surprising, since the ambient magnetic field lines in the 
low corona are almost potential, even in the MHD solution. 
The main differences between the MHD solution and the 
potential field appear in the field lines that are stretched by 
the solar wind in the high corona and in the field lines of the 
superimposed fiux rope, which are twisted. Unlike the po- 
tential field method, the time-dependent, MHD simulation 
enables us to follow the temporal evolution of the magnetic 
field and its dynamic response to the CME expansion. 

Figure 3 shows snapshots of the evolution of the field 
topology from the start of the simulation (row 1) to one 
(row 2) and four (row 3) hours later. The left column shows 
the photosphere colored with contours of Br intensity along 
with selected magnetic field lines (which are not the same 
in all frames). Magenta marks field lines in the core CME 
flux rope, with both footpoints located at the source region. 
Cyan marks overlying field lines with one footpoint rooted 
in the source region, and yellow marks field lines of the over- 
lying helmet streamers and open field lines with footpoints 
located far from the source region. Row 2 shows that after 
one hour, some legs of the three-dimensional CME (repre- 
sented by the highly twisted field lines in cyan, magenta or 
yellow) have expanded beyond the original active region so 
that there are more cyan field lines and fewer magenta ones. 
After four hours, row 3 shows that many field lines of the 
CME are located far from the source active region, though 
we emphasize that some still remain rooted at the original 
source AR. While the CME dynamic expansion pushes the 
surrounding field lines to the sides, the only way the foot- 
point of the CME can migrate away from the source AR 
is through reconncction with surrounding field lines as the 
CME is expanding. By the time the CME reaches the height 
of the large helmet streamers and the open fiux that have 
opposite polarity to that of the CME, the CME footpoint 
expansion stops and only the traditional volumetric expan- 
sion of the CME continues. The overall field evolution at 
this point in the simulation is roughly consistent with the 
idealized scenario suggested by the breakout model, proba- 
bly due to the similar field topologies of the idealized and 
real cases. 

The right column of Figure 3 shows a sphere at height 
r = l.O67?0 colored with the sign of Br for the initial state 
of the simulation, where positive is yellow and negative is 
blue. Overlain on this initial polarity pattern is the time- 
evolving polarity inversion line in red. Rows 2 and 3 show 
that the inversion line changes as the CME propagates and 
interacts with the surrounding field. This change in the dis- 
tribution of positive-negative flux reflects the exchange of 
closed flux between the CME and the surrounding field. 

In order to validate the simulated CME flux rope as it 
propagates into interplanetary space, we extract the sim- 
ulated plasma parameters across the flux rope (upstream 



to downstream) from the frame of t = lOh after the erup- 
tion onset. We cannot directly compare our simulation re- 
sult with 1 AU in-situ data, since the simulated CME only 
reaches ~ 23 — 2ARq , and the coordinate systems of the two 
data sets are different (Heliocentric Caxrington system in 
the simulation vs. Geocentric system for in-situ data) . How- 
ever, we can compare the fiux rope orientation and chirality 
through the local Bz component and assume these should 
be conserved from the simulated location of the CME to 1 
AU. 

The first three panels of Figure 4 show a line extrac- 
tion of S, Br, and Bz along constant latitude and longi- 
tude from upstream (un-shocked plasma) to downstream 
(shocked plasma) in the simulation frame of reference. The 
fourth panel shows the GSE Bz component observed by 
the Wind spacecraft during May 14-16 (data taken from 
http://cdaweb.gsfc.nasa.gov/). It can be seen that the field 
rotation, indicated by the change of sign in Bz, is present in 
the simulation as in the observations (compare panels three 
and four) . However, we note that the magnitude of the sim- 
ulated Bz is greater than in the observations. This may be 
attributed to the simulation data being extracted closer to 
the Sun than the 1 AU observations. Here we do not attempt 
to match the result of the model with this observation, which 
has been taken at different location in the heliosphere and 
in different coordinate system, only to demonstrate that the 
field rotation seen in the data also appears in the simulation. 

A visualization of the field topology (also showing the 
line of extraction) is displayed in Figure 5. The left panel 
shows the photosphere as the inner sphere colored by radial 
field strength. The approximate CME front is shown by the 
gray shading, which outlines an iso-surface of density ratio 
equal to 4 between the t = lOh frame and the initial state. 
Selected field lines illustrate the field topology. Blue and 
red indicate negative and positive values of Bz , respectively, 
and the cartoon of a hand indicates the left-handed chirality 
of the flux rope, where the thumb points in the direction of 
the axial field and the fingers curl in the direction of the 
coils. These are in agreement with the line plots of Bz in 
Figure 4, representing a transition from a downward field 
(negative Bz) at the front of the CME to an upward field 
(positive Bz) at the back. The signature is a smooth ro- 
tation of the field direction, indicating the passage across 
the core of a flux rope. This is one of the most distinctive 
signatures of a magnetic cloud [Burlaga et al., 1981]. The 
CME topology in the left panel of Figure 5 is also consistent 
with the observed interplanetary magnetic cloud shown in 
the white inset in the lower right corner (from Crooker et al. 
[2008]). The right panel of Figure 5 shows a view looking 
down on the ecliptic plane, which is colored with contours of 
the solar wind speed. Selected field lines are shown in black, 
and the blue arrow indicates the line of extraction used for 
the plots in Figure 4. 

Figure 6 shows the orientation of the fiux rope at differ- 
ent stages of the simulation, where the display is similar to 
the left panel of Figure 5. It can be seen that the rotation 
is well established about 6 hours after the eruption onset 
and that, after this, the CME maintains its orientation (it 
can still rotate around the z-axis) . This result is also consis- 
tent with the overall behavior of the flux rope rotation de- 
scribed by Lynch et al. [2009], who used the breakout model 
in a topologically idealized numerical simulation to study the 
flux rope rotation. In their model, a rotation of ~ 40 — 50 
degrees is established during the first ~ 20 minutes of the 
flare onset (co-temporal with an increase in the magnitude 
of the poloidal magnetic field). They ascribe the rotation 
primarily to the Lorentz force and overall torque from the 



COHEN ET AL.: INTERCHANGE RECONNECTION IN THE CORONA 



X- 5 



sheared core magnetic field and related tension forces and 
conclude that the kink instability [e.g. Kliem et al, 2004], 
[see also Green et al., 2007] is not responsible. We cannot 
study the possible effects of the kink instability in our sim- 
ulation because the Titov-Demoulin flux rope is idealized 
and no kink-unstable twist is introduced. Our simulation 
results seem to imply that the flux rope rotation is strongly 
influenced by interaction with the surrounding field, since 
the rotation is established at the early stage of the CME 
propagation, between 04:00 and 06:00 li (Figure 6). 

In Figure 7, wc display selected field lines and show their 
evolution in time in order to illustrate the stepping rccon- 
nection between the CME and adjacent loops, as described 
in Section 1. Unlike the field lines in Figure 3, here the 
same field lines are displayed in the different time frames. 
The sphere represents the photosphere colored with the den- 
sity base differences as in the right panel of Figure 2. In the 
left and middle panels of Figure 7, the red line indicates 
the flux rope and has both footpoints rooted in the source 
AR. The cyan field line represents a neighboring field line 
with only one footpoint rooted in the source AR. The CME 
initially expands with its flux rope footpoints (red) located 
in the active region (middle panel). Later on (right panel), 
reconnection has occurred between the expanding flux rope 
(red), and the neighboring field line (cyan). This reconnec- 
tion transfers part of the original twisted flux rope field to 
the neighboring loop. The expanded flux rope now has only 
one footpoint in the source active region. The other has 
been "stepped out" and is now rooted in the quiet Sun at a 
considerable distance from the active region. The results dis- 
played in Figure 7 give an example of how one field line in a 
CME expands laterally as a result of magnetic reconnection 
between the core flux rope and surrounding magnetic field. 
This mechanism migrates many of the footpoints of the field 
lines associated with the CME away from the source AR and 
thus facilitates a global lateral expansion of the CME in the 
low corona [van Dnel-Gesztelyi et al., 2008]. 

4. Discussion 

The results of our numerical simulation show that in the 
low corona the CME interacts with the surrounding field, 
which is composed of magnetic loops with different sizes, 
orientations, and mixed polarity. The path of CME expan- 
sion in this region is determined by the easiest way the CME 
can get through the surrounding field structure. This path 
depends on where the magnetic tension is smallest (i.e., in 
larger overlying loops) and where the field topology allows 
the CME to reconnect with surrounding field lines. The lat- 
ter leads to a lateral expansion of the CME via migration 
of the CME footpoints. In addition, the simulation results 
show that the left-handed flux rope rotates counterclockwise 
by 90 degrees by the time the CME has expanded into in- 
terplanetary space so that the original ~ N-S orientation of 
the flux rope (see Figure 1), is now aligned ~ E-W. This is 
consistent with the orientation of the magnetic cloud in the 
ecliptic plane deduced from observations [Webb et al., 2000] 
and from modeUng [Attnll et al, 2006). Both the CME lat- 
eral expansion and rotation wane when the CME starts to 
interact with the more organized magnetic flux of the large- 
scale helmet streamers and open flux. 

It is reasonable to believe that since both the CME ro- 
tation and expansion in our simulation arc the result of the 
interaction between the expanding flux rope and the ambi- 
ent flux via magnetic reconnection, the sense of the rotation 
is related to the chirality of the flux rope as well as the ori- 
entation of the ambient field, as suggested by Lynch et al. 
[2005]. We further suggest that the final orientation of the 
flux rope represents the "relaxed" state of the reconnection 



process between the core flux rope and the surrounding coro- 
nal field. 

Overall, our simulation produces results that arc consis- 
tent with the breakout model, primarily because the topol- 
ogy of the ambient field and the core flux rope for the May 
12, 1997, event is similar to the idealized field topology in 
the breakout model. One should keep in mind, however, 
that this might not be the case during solax maximum, when 
both the ambient field and the orientation of the core flux 
rope are unlikely to fulfill such idealized conditions. In ad- 
dition, we stress that in the breakout model the field topol- 
ogy drives the eruption, while our simulation captures the 
post-eruption interaction between the CME and the ambient 
field. 

The result of the simulation suggests a scenario which en- 
ables us to develop a more complete picture of the May 12, 
1997, eruption. Panels (la) and (lb) of Figure 8 illustrate 
the scenario proposed earlier by Crooker and Webb [2006]; 
Attrill et al. [2006], and Crooker et al. [2008], as well as 
Owens et al. [2007]. The CME, represented by the twisted 
line in panel (la), expands and eventually undergoes inter- 
change reconnection (ICX) with an open field line of the 
north polar coronal hole. As a result, the magnetic field 
in the CME that is rooted in the south dimming region 
is essentially an open field line (panel lb), consistent with 
the open fields deduced from electron observations in the 
magnetic cloud at 1 AU. The accumulation of progressively 
larger closed loops between the northern dimming region 
and the polar coronal hole has the effect of contracting the 
northern dimming region, which is well underway by 08:26 
UT (see Figure 4 in Attnll et al. [2006]). The south dimming 
region remains near its maximum spatial extent for a longer 
time than the north dimming region, but by 15:00 UT, it 
has started to contract, as well. This contraction was not 
clearly addressed by Crooker and Webb [2006] and Attrill 
et al. [2006]. Why should the south dimming region should 
show a systematic contraction when the field lines rooted in 
that region should remain open (Figure 8, panel lb)? 

Our simulation suggests a slightly different scenario, that 
provides an answer to that question. In this scenario, illus- 
trated in panels 2a and 2b of Figure 8, the core flux rope first 
reconnects with greatly extended overlying closed field lines. 
As a result, the legs of the twisted closed loops that consti- 
tute the CME are displaced so that only the south dimming 
region remains connected to the greatly expanded twisted 
CME loops (panel 2b, Figure 8). We note that a flux rope 
is made up of many complex twisted magnetic field lines. 
The outermost field lines reconnect with the surrounding 
magnetic field, but the innermost core field lines may not 
experience such reconnections. As a result, some part of the 
original flux rope may remain rooted in the source region 
throughout the eruption (not shown in Figure 8). Indeed, 
such a picture is consistent with the deepest parts of coro- 
nal dimming regions remaining dimmed for several days [e.g. 
Attrill et al, 2008], long after the main bulk of a CME has 
left the Sun. 

The contraction of the south dimming region in this May 
12, 1997, event can now be explained as a result of suc- 
cessively larger closed loops being formed duo to ongoing 
breakout-type reconnections, so contracting the south dim- 
ming in a similar manner to the north one, albeit at a later 
time due to the much larger scale of those closed loops. Dur- 
ing this process, one of the legs of the twisted field making 
up the CME is displaced out of the source AR. To our knowl- 
edge, this is the first time that a CME has been shown to 
comprise one of the side lobes created as a result of reconnec- 
tion with the overlying field. Usually the central flux system 
in the breakout model erupts as the restraining overlying 
field is removed to the side lobes. Here the situation differs 
markedly from the textbook breakout scenario in several 
ways: (i) the reconnection occurs much closer to one foot- 
point of the overlying closed field than the other, introducing 
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a marked asymmetry; and (ii) the overlying elosed field is 
greatly extended so that upon reconnection, relatively long 
field lines constitute the side lobe and become dragged out 
into interplanetary space as the eruption proceeds. These 
two factors mean that the north dimming region will start 
to contract, whilst the south region will remain extended for 
a longer time, as observed. 

What the new scenario in Figure 8 may struggle to ex- 
plain is the unidirectional streaming of suprathermal elec- 
trons in the CME at 1 AU, indicating field lines connected 
to the Sun at only one end. Although the determination of 
magnetic connectivity from suprathermal electron data can 
be hindered by significant uncertainties [e.g. Riley et al., 
2004], in this case the signature of open fields throughout 
most of CME is quite clear. In contrast, our simulation re- 
sults show that the extended twisted field identified as the 
CME is connected to the Sun at both ends, which topo- 
logically would be expected to produce counterstreaming 
electrons at 1 AU. Two possible explanations for this dis- 
agreement are: 1) Although not apparent in the simulation, 
the reconnection with overlying loops pictured in panel lb of 
Figure 8 may proceed until the supply of closed loops gives 
way to the neighboring open polar field, in which case the 
scenario essentially reverts to the original one where inter- 
change reconnection in the corona opens the loops (panel 
2c, Figure 8); 2) interchange reconnection opens the loops 
associated with the CME, but takes place gradually, in the 
solar wind, outside the simulation box, in transit to 1 AU. 

The first explanation suggests that the helmet streamer 
(forming the initial overlying closed magnetic field) is tem- 
porarily removed by the breakout process, and then later 
reforms through ICX with the open fiux of the polar region. 
Even though we believe that the ICX scenario represented 
by panel 2c in Figure 8 must occur at some point, it is hard 
to capture this interaction in the current simulation (in con- 
trast to the interaction with the closed fiux). This may be 
due to the fact that the ICX interaction should occur closer 
to the outer boundary, where every field line that leaves the 
simulation domain is considered to be "open". We believe 
that the interaction of the CME flux with the open flux 
of the Sun can be better studied through idealized simula- 
tions, in which the isolation of such reconnection events can 
be better resolved, rather than in this complicated, realistic 
simulation. In addition, we favor the first explanation over 
the second because reconnection in the solar wind does not 
act to balance the fiux increase owing to the introduction 
of the CME loops [e.g. Owens and Crooker, 2006], nor does 
it transport open flux on the Sun in the course of the solar 
cycle [e.g. Owens et al., 2007]. At this stage, however, we 
leave the issue as an open question. 

5. Summary and Conclusions 

Using numerical simulation, we study the evolution of 
a CME and its interaction with the coronal magnetic field 
based on the May 12, 1997, CME event. Our simulation 
provides the following results: 

• The simulated density corresponding to a cut across 
the CME at 1.1_R0 is comparable to the coronal wave bright 
front that appears in the observed EIT images. This is sim- 
ilar to results obtained for another event studied by Cohen 
et al. [2009]. 

• During the initial stage of CME expansion, the core flux 
rope reconnects with the surrounding (neighboring) field, 
which facilitates lateral expansion of the CME footprint in 
the low corona. Evidence for this is the displacement of 
the initial CME footpoints away from the source AR, via 
"stepping reconnection" . 

• The CME then reconnects with the oppositely orien- 
tated overlying magnetic field in the manner of the breaJf- 
out model. This is due to the global field topology, which 
is essentially the same as the initial state of the breakout 
model. 



• A significant component of the CME that expands into 
interplanetary space comprises one of the side lobes created 
mainly as a result of reconnection with the overlying field. 

• As the flux rope expands, it rotates counterclockwise 
by 90 degrees owing not to the kink instability but to recon- 
nection with the surrounding fields. 

• The lateral expansion as well as the rotation of the flux 
rope due to interaction with the overlying magnetic field 
continue as long as the reconnections occur. They stop when 
the CME reaches the helmet streamers and the open flux, 
where the field of the CME matches the globally-organized 
field. 

• The orientation and left-handed chirality of the simu- 
lated flux rope are in agreement with the interplanetary flux 
rope observed at 1 AU. 

• The simulation shows no direct evidence of the antici- 
pated interchange reconnection that would account for the 
electron observations indicating open fields in the CME at 
1 AU. 

From these results we conclude that reconnection between 
closed loops may play a major role in transforming a small- 
scale flux rope in an active region to a large-scale flux rope 
with widely-spaced footpoints and a new orientation before 
leaving the Sun as a CME. Therefore, this complex interac- 
tion should be talcen into consideration in any attempt to 
predict the CME geo-efFectiveness. We further conclude that 
simulations of CMEs that erupt into realistic background 
fields constructed from magnetograms afford unprecedented 
views of the complicated three-dimensional development of 
magnetic structure. 

Acknowledgments. 

We thank for two referees for their useful comments and 
suggestions. OC and NUC are supported by SHINE through 
NSF grants ATM-0823592 and ATM-0553397. NAS was sup- 
ported through the NASA LWS EMMREM project, grant 
no. NNX07AC14G. GDRA is supported by NASA grant 
NNX09AB11G-R. Simulation results were obtained using the 
Space Weather Modelling Framework, developed by the Center 
for Space Environment Modelling, at the University of Michigan 
with funding support from NASA ESS, NASA ESTO-CT, NSF 
KDI, and DoD MURI. SOHO is a project of international coop- 
eration between ESA and NASA. 



References 

Altschulcr, M. D., and G. Ncwkirk, Magnetic Fields and the 
Structure of the Solar Corona. 1: Methods of Calculating Coro- 
nal Fields, Sol. Phys., 9, 131-149, doi;10.1007/BF00145734, 
1969. 

Amari, T., J. F. Luciani, Z. Mikic, and J. Linker, A Twisted 
Flux Rope Model for Coronal Mass Ejections and Two- Ribbon 
Flares, Astrophys. J., , 529, L49-L52, doi:10.1086/312444, 
2000. 

Antiochos, S. K., C. R. DeVore, and J. A. Klimchuk, A Model for 
Solar Coronal Mass Ejections, Astrophys. J., , 510, 485-493, 
doi: 10. 1086/306563, 1999. 

Arge, C. N., and V. J. Pizzo, Improvement in the Prediction 
of Solar Wind Conditions Using Near-Real time Solar Mag- 
netic Field Updates, Journal of Geophysical Research (Space 
Physics), 105, 10,465-10,479, doi: 10. 1029/1999 JA900262, 
2000. 

Attrill, C, M. S. Nakwacki, L. K. Harra, L. van Driel-Gesztelyi, 
C. H. Mandrini, S. Dasso, and J. Wang, Using the Evolution 
of Coronal Dimming Regions to Probe the Global Magnetic 
Field Topology, Sol. Phys., 238, 117-139, doi:10.1007/sll207- 
006-0167-5, 2006. 



COHEN ET AL.: INTERCHANGE RECONNECTION IN THE CORONA 



X- 7 



Attrill, G. D. R., L. K. Harra, L. van Dricl-Gcsztclyi, and 
P. Dcmoulin, Coronal "Wave" : Magnetic Footprint of a Coro- 
nal Mass Ejection?, Astrophys. J., , 656, L101-L104, doi: 
10.1086/512854, 2007. 

Attrill, G. D. R., L. van Dricl-Gesztelyi, P. Dcmoulin, A. N. 
Zliukov, K. Steed, L. K. Harra, C. H. Mandrini, and J. Linker, 
The Recovery of CME-Rclated Dimmings and the ICME's En- 
during Magnetic Connection to the Sun, Sol. Phys., 252, 349— 
372, doi;10.1007/sll207-008-9255-z, 2008. 

Attrill, G. D. R., A. .J. EngelL M. .1. Wills-Davey, P. Grigis, 
and P. Testa, Hinode/XRT and STEREO Observations of 
a Diffuse Coronal "Wave" -Coronal Mass Ejection-Dimming 
Event, Astrophys. J., , 704, 1296-1308, doi: 10. 1088/0004- 
637X/704/2/1296, 2009. 

Biesecker, D. A., D. C. Myers, B. J. Thompson, D. M. Hammer, 
and A. Vourlidas, Solar Phenomena Associated with "EIT 
Waves", Astrophys. J., , 569, 1009-1015, doi: 10. 1086/339402, 
2002. 

Boyd, T. J. M., and J. J. Sanderson, The Physics of Plasmas, 
2003. 

Burlaga. L., E. Sittler, F. Mariani, and R. Schwenn, J. Geophys. 
Res., , 86, 6673, 1981. 

Cohen, O., I. V. Sokolo, I. I. Roussev, N. Lugaz, W. B. Manch- 
ester, T. I. Gombosi, and C. N. Arge, Validation of a global 
3D heliospheric model with observations for the May 12, 1997 
CME event, Journal of Atmospheric and Solar- Terrestrial 
Physics, 70, 583-592, doi:10.1016/j.jastp.2007.08.065, 2008a. 

Cohen, C, I. V. Sokolov, 1. 1. Roussev, and T. I. Gombosi, Valida- 
tion of a synoptic solar wind model, Journal of Geophysical Re- 
search (Space Physics), 113, 3104, doi:10.1029/2007JA012797, 
2008b. 

Cohen, C, G. D. R. Attrill, W. B. Manchester, and M. J. 
Wills-Davey, Numerical Simulation of an EUV Coronal Wave 
Based on the 2009 February 13 CME Event Observed by 
STEREO, Astrophys. J., , 705, 587-602, doi:10.1088/0004- 

637X/705/1/587, 2009. 
Cohen, O., et al., A Semi-Empirical Magnetohydrodynamical 
Model of the Solar Wind, Astrophys. J., , 645, L163-L166, 
2007. 

Crooker, N. U., and D. F. Webb, Remote sensing of the solar site 
of interchange reconncction associated with the May 1997 mag- 
netic cloud. Journal of Geophysical Research (Space Physics), 
111, 8108, doi:10.1029/2006JA011649, 2006. 

Crooker, N. U., J. T. Gosling, and S. W. Kahler, Reducing helio- 
spheric magnetic flux from coronal mass ejections without dis- 
connection. Journal of Geophysical Research (Space Physics), 
107, 1028, doi:10.1029/2001JA000236, 2002. 

Crooker, N. U., S. W. Kahler, J. T. Gosling, and R. P. Lepping, 
Evidence in magnetic clouds for systematic open flux transport 
on the Sun, Journal of Geophysical Research ( Space Physics ), 
113, 12,107, doi:10.1029/2008JA013628, 2008. 

Delannee, C, The role of small versus large scale magnetic topolo- 
gies in global waves, Astron. Astrophys., , 4^5, 571—575, doi: 
10.1051/0004-6361:200809451, 2009. 

Dere, K. P., et al., EIT and LASCO Observations of the Initi- 
ation of a Coronal Mass Ejection, Sol. Phys., 175, 601-612, 
doi: 10. 1023/A: 1004907307376, 1997. 

Downs, C, I. I. Roussev, B. van der Hoist, N. Lugaz, I. V. 
Sokolov, and T. I. Gombosi, Toward a Realistic Thermo- 
dynamic Magnetohydrodynamic Model of the Global Solar 
Corona, Astrophys. J., , 712, 1219-1231, doi: 10. 1088/0004- 
637X/712/2/1219, 2010. 

Fan, Y., and S. E. Gibson, Onset of Coronal Mass Ejections Due 
to Loss of Confinement of Coronal Flux Ropes, Astrophys. J., 
, 668, 1232-1245, doi:10. 1086/521335, 2007. 

Forbes, T. G., and E. R. Priest, Photospheric Magnetic Field 
Evolution and Eruptive Flares, Astrophys. J., , 44^^ 377, doi: 
10.1086/175797, 1995. 

Forbes, T. G., et al., CME Theory and Models, Space Science 
Reviews, 123, 251-302, doi:10.l607/sll214-006-9019-8, 2006. 

Green, L. M., B. Kliem, T. Torok, L. van Dricl-Gesztelyi, 
and G. D. R. Attrill, Transient Coronal Sigmoids and Ro- 
tating Erupting Flux Ropes, Sol. Phys., 246, 365—391, doi: 
10.1007/sll207-007-9061-z, 2007. 

Kliem, B., V. S. Titov, and T. Torok, Formation of current sheets 
and sigmoidal structure by the kink instability of a magnetic 
loop, Astron. Astrophys., , 413, L23-L26, doi: 10. 1051/0004- 
6361:20031690, 2004. 



Landi, E., U. Feldman, and K. P. Dere, CHIANTI-An Atomic 
Database for Emission Lines. V. Comparison with an Isother- 
mal Spectrum Observed with SUMER, Astrophys. J. (Supp.), 
, 139, 281-296, doi:10. 1086/337949, 2002. 

Linker, J. A., Z. Mikic, R. Lionello, P. Riley, T. Amari, and 
D. Odstrcil, Flux cancellation and coronal mass ejections. 
Physics of Plasmas, 10, 1971-1978, doi:10.1063/l. 1563668, 
2003. 

Lugaz, N., W. B. Manchester, IV, I. I. Roussev, G. Toth, and T. I. 
Gombosi, Numerical Investigation of the Homologous Coronal 
Mass Ejection Events from Active Region 9236, Astrophys. ./., 
, 659, 788-800, doi: 10. 1086/512005, 2007. 

Lynch, B. J., J. R. Gruesbeck, T. H. Zurbuchen, and S. K. An- 
tiochos. Solar cycle-dependent helicity transport by magnetic 
clouds, Journal of Geophysical Research (Space Physics), 110, 
8107, doi:10.1029/2005JA011137, 2005. 

Lynch, B. J., S. K. Antiochos, C. R. DeVore, J. G. Luhmann, 
and T. H. Zurbuchen, Topological Evolution of a Fast Mag- 
netic Breakout CME in Three Dimensions, Astrophys. J., , 
683, 1192-1206, doi:10.1086/589738, 2008. 

Lynch, B. J., S. K. Antiochos, Y. Li, J. G. Luhmann, and 
C. R. DeVore, Rotation of Coronal Mass Ejections during 
Eruption, Astrophys. J., , 697, 1918-1927, doi: 10. 1088/0004- 
637X/697/2/1918, 2009. 

Manchester, W. B., T. I. Gombosi, I. Roussev, A. Ridley, D. L. 
De Zeeuw, I. V. Sokolov, K. G. Powell, and G. Toth, Mod- 
eling a spax;e weather event from the Sun to the Earth: 
CME generation and interplanetary propagation. Journal 
of Geophysical Research (Space Physics), 109, 2107, doi: 
10.1029/2003JA010150, 2004. 

Manchester, W. B., IV, A. Vourlidas, G. Toth, N. Lugaz, I. 1. 
Roussev, I. V. Sokolov, T. I. Gombosi, D. L. De Zeeuw, and 
M. Opher, Three-dimensional MHD Simulation of the 2003 
October 28 Coronal Mass Ejection: Comparison with LASCO 
Coronagraph Observations, Astrophys. J., , 684, 1448-1460, 
doi:10.1086/590231, 2008. 

Manoharan, P. K., L. van Dricl-Gesztelyi, M. Pick, and P. De- 
moulin, Evidence for Large-Scale Solar Magnetic Reconnection 
from Radio and X-Ray Measurements, Astrophys. J., , .^SS, 
L73-t-, doi:10.1086/310221, 1996. 

Moses, D., et al., EIT Observations of the Extreme Ultraviolet 
Sun, Sol. Phys., 175, 571-599, doi:10.1023/A:1004902913117, 
1997. 

Odstrcil, D., P. Riley, and X. P. Zhao, Numerical simulation 
of the 12 May 1997 interplanetary CME event. Journal of 
Geophysical Research (Space Physics), 109, 2116 — h, doi: 
10.1029/2003JA010135, 2004. 

Odstrcil, D., V. J. Pizzo, and C. N. Arge, Propagation of the 
12 May 1997 interplanetary coronal mass ejection in evolving 
solar wind structures. Journal of Geophysical Research (Space 
Physics), 110, 2106, doi:10.1029/2004JA010745, 2005. 

Ofman, L., Three-dimensional MHD Model of Wave Activity in a 
Coronal Active Region, Astrophys. J., , 655, 1134—1141, doi: 
10.1086/510012, 2007. 

Ofman, L., and B. J. Thompson, Interaction of EIT Waves with 
Coronal Active Regions, Astrophys. ./. , , 574, 440-452, doi: 
10.1086/340924, 2002. 

Owens, M. J., and N. U. Crooker, Coronal mass ejections 
and magnetic flux buildup in the lieliospliere, Journal of 
Geophysical Research (Space Physics), 111, 10,104, doi: 
10.1029/2006JA011641, 2006. 

Owens, M. J., N. A. Schwadron, N. U. Crooker, W. J. Hughes, 
and H. E. Spence, Role of coronal mass ejections in the he- 
liospheric Hale cycle, Geophys. Res. Lett., , 34 , 6104, doi: 
10.1029/2006GL028795, 2007. 

Pick, M., et al.. Joint Nancay Radioheliograph and LASCO Ob- 
servations of Coronal Mass Ejections - II. The 9 July 1996 
Event, Sol. Phys., 181, 455-468, 1998. 

Podladcliikova, O., and D. Bcrghmans, Automated Detection 
Of Eit Waves And Dimmings, Sol. Phys., 228, 265-284, doi: 
10.1007/sll207-005-5373-z, 2005. 

Pohjolainen, S., et al., On-the-Disk Development of the Halo 
Coronal Mass Ejection on 1998 May 2, Astrophys. J., , 556, 
421-431, doi:10.1086/321577, 2001. 

Powell, K. G., P. L. Roe, T. J. Linde, T. I. Gombosi, and D. L. 
de Zeeuw, A Solution- Adaptive Upwind Scheme for Ideal Mag- 
netohydrodynamics, Journal of Computational Physics, 154, 
284-309, doi:10.1006/jcph.l999.6299, 1999. 



X- 8 



COHEN ET AL.: INTERCHANGE RECONNECTION IN THE CORONA 



Riley, P., J. T. Gosling, and N. U. Crooker, Ulysses Observa- 
tions of the Magnetic Connectivity between Coronal Mass 
Ejections and the Sun, Astrophys, J,, , 608, 1100-1105, doi: 
10.1086/420811, 2004. 

Riley, P., R. Lionello, Z. Mikic, and J. Linker, Using Global Sim- 
ulations to Relate the Three-Part Structure of Coronal Mass 
Ejections to In Situ Signatures, Astrophys. J., , 672, 1221- 
1227, doi:10.1086/523893, 2008. 

Roussev, I. I., T. G. Forbes, T. I. Gombosi, I. V. Sokolov, D. L. 
DeZeeuw, and J. Birn, A Three-Dimensional Flux Rope Model 
for Coronal Mass Ejections Based on a Loss of Equilibrium, 
Astrophys. J., , 588, L45-L48, 2003a. 

Roussev, I. I., T. I. Gombosi, I. V. Sokolov, M. Velli, W. Manch- 
ester, D. L. DeZccuw, P. Liewer, G. Toth, and J. Luhmann, A 
Three-Dimensional Model of the Solar Wind Incorporating So- 
lar Magnctogram Observations, Astrophys. J., , 595, L57-L61, 
2003b. 

Schmidt, J. M., and L. Ofman, Global Simulation of an Ex- 
treme Ultraviolet Imaging Telescope Wave, Astrophys. -J., , 
713, 1008-1015, doi;10.1088/0004-637X/713/2/1008, 2010. 

Shen, P., X. Feng, S. T. Wu, and C. Xiang, Three-dimensional 
MHD simulation of CMEs in three-dimensional background so- 
lar wind with the self-consistent structure on the source surface 
as input: Numerical simulation of the January 1997 Sun-Earth 
connection event. Journal of Geophysical Research (Space 
Physics), 112, 6109— h, doi:10.1029/2006JA012164, 2007. 

Sokolov, I. v., K. G. Powell, O. Cohen, and T. I. Gombosi, Com- 
putational MagnetoHydrodynamics, Based on Solution of the 
Well- Posed Riemann Problem, in Numerical Modeling of Space 
Plasma Flows, Astronomical Society of the Pacific Confer- 
ence Series, vol. 385, edited by N. V. Pogorelov, E. Audit, & 
G. P. Zank, pp. 291-+, 2008. 

Thompson, B. J., S. P. Plunkctt, J. B. Gurman, J. S. Newmark, 
O. C. St. Cyr, and D. J. Michels, SOHO/EIT observations 



of an Earth-directed coronal mass ejection on May 12, 1997, 
Geophys. Res. Lett, , 25, 2465-2468, doi:10.1029/98GL50429, 
1998. 

Titov, V. S., and P. Dcmoulin, Basic Topology of Twisted Mag- 
netic Configurations in Solar Flares, Astron. Astrophys. , , 351, 
707-720, 1999. 

Toth, G., ct al., Space Weather Modeling Framework; A New 
Tool for the Space Science Community, J. Geophys. Res., 110, 
12,226-12,237, doi;10.1029/2005.IA011126, 2005. 

van dcr Hoist, B., W. Manchester, I. V. Sokolov, G. Toth, T. I. 
Gombosi, D. DeZecuw, and O. Cohen, Breakout Coronal Mass 
Ejection or Streamer Blowout; The Bugle Effect, Astrophys. 
J., , 693, 1178-1187, doi;10.1088/0004-637X/693/2/1178, 
2009. 

van Driel-Gesztclyi, L., C. P. Goff, P. Dcmoulin, ,1. L. Cul- 
hanc, S. A. Matthews, L. K. Harra, C. H. Mandrini, 
K. Klein, and H. Kurokawa, Multi-scale reconnections in a 
complex CME, Advances in Space Research, 42, 858—865, doi: 
10. 1016/j.asr. 2007.04.065, 2008. 

Wang, Y.-M., and N. R. Sheeley, Solar Wind Speed and Coronal 
Flux- Tube Expansion, Astrophys. J., , 355, 726-732, 1990. 

Webb, D. F., R. P. Lepping, L. F. Burlaga, C. E. DeForest, D. E. 
Larson, S. F. Martin, S. P. Plunkett, and D. M. Rust, The ori- 
gin and development of the May 1997 magnetic cloud, J. Geo- 
phys. Res., , 105, 27,251-27,260, doi;10.1029/2000JA000021, 
2000. 

Wu, C, C. D. Fry, S. T. Wu, M. Dryer, and K. Liou, Three- 
dimensional global simulation of interplanetary coronal mass 
ejection propagation from the Sun to the heliospherc: Solar 
event of 12 May 1997, -Journal of Geophysical Research (Space 
Physics), 112, 9104— 1-, doi:10.1029/2006JA012211, 2007. 



COHEN ET AL.: INTERCHANGE RECONNECTION IN THE CORONA X - 9 




Figure 1. The left panel shows the large-scale structure 
of the steady-state corona. Selected field lines are shown 
in blue (closed field lines) and red (open field lines). 
Background color contours represent the solar wind ra- 
dial speed (yellow- fast, blue-slow). The right panel shows 
the vicinity of the active region at the initial state of the 
simulation, where its approximated location on the Sun 
is marked by the black square. Color contours represent 
the magnitude of the radial field on a sphere at a height 
of 1.06-Ro, red streamlines represent three-dimensional 
magnetic field lines of the superimposed flux rope, and 
solid white lines mark the grid structure around the flux 
rope. The flux rope itself is represented by an iso surface 
of mass density, p = 2 • 10~^* g cm~^, which is greater 
than the surrounding density at the same height. 
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Figure 2. Comparison of observed and simulated coro- 
nal dimmings and the coronal wave for the early stage of 
the simulation: 10, 30, 40, and 60 minutes after the erup- 
tion onset (top to bottom). Left column shows SOHO 
EIT 195A base-difference images, second column shows 
the corresponding synthetic EIT 195^ base-difference im- 
ages produced by the simulation. The third panel shows 
base-difference images of the simulated mass density on a 
sphere at height of r = 1.1_Rq. The right column shows a 
similar display to the third panel, but with the addition of 
selected magnetic field lines of the core flux rope marked 
in red, selected field lines with one footpoint rooted in the 
original source region in blue, and selected overlying field 
lines with both footpoints located far from the source re- 
gion in yellow. The white circle in the first three columns 
indicates the approximate location of the leading edge of 
the coronal wave disturbance. 
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Figure 3. Evolution of the field topology at later time of 
the simulation than shown in Figure 2. Left panel shows 
the photosphere colored with contours of Br intensity, 
along with selected magnetic field lines. Magenta field 
lines correspond to the core CME flux rope, with both 
footpoints rooted in the source region. Cyan lines mark 
overlying field with one footpoint located at the source 
region, and yellow lines mark field lines of the overlying 
helmet streamers and open field lines with footpoints lo- 
cated far from the source region. Right panel shows a 
sphere at a height of r = 1.06-Rq colored with the sign of 
Br (positive-yellow, negative-blue) for the initial state of 
the simulation. The red line marks the inversion line for 
the particular time frame. 
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Figure 4. First three panels sliow upstream to down- 
stream line extraction of B, Br, and B^, respectively, 
in the simulation frame of reference. The fourth panel 
shows the GSE B^ observed by the Wind spacecraft be- 
tween May 14-16. 
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Figure 5. The field topology at t = lOh. Left panel 
shows the photosphere as the inner sphere colored with 
the radial field line contours. The approximate CME 
front is shown by the white shade, which indicates an iso- 
surface of density ratio of 4 between the t = lOh frame 
and the initial state. Selected field lines show the front 
of the CME and the field topology. Blue color indicates 
negative value of , red color indicates positive value of 
Bz, and the left hand cartoon represents the flux rope 
chirality. Inset shows the observed interplanetary mag- 
netic cloud chirality and orientation (from Crooker et al. 
[2002] . Right panel shows a view on the equatorial plane 
colored with contours of the solar wind speed together 
with selected three-dimensional field lines. The blue ar- 
row marks the upstream-downstream line of extraction 
shown in Figure 4 
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Figure 6. The orientation of tlie flux rope at different 
stages of the simulation. Display is similar to the left 
panel of Figure 5. Blue (red) indicates negative (positive) 
Bz. The orientation of the flux rope axis is initially ~N-S 
(cf. Figure 1). By 00:40 h, the orientation of the axis is 
NW-SE. In the frame at 04:00 h, the orientation is hard 
to establish, but by 06:00 h, the axis of the flux rope lies 
~E-W. The last two frames at 12:00 h and 14:00 h are 
taken when most of the CME has left the simulation box. 




Figure 7. Displacement of the CME footpoints via step- 
ping reconnection with neighboring magnetic field. The 
figure shows selected field lines and their evolution in 
time. Density base differences at a height of r = I.IRq 
are shown in grey (as in Figure 2). In the left and middle 
panels, the red line indicates the flux rope and has both 
footpoints rooted in the source AR. The cyan field line 
has only one footpoint rooted in the source AR, and rep- 
resents a neighboring field line. In all frames, twisted field 
lines are associated with the CME. The CME initially ex- 
pands with its flux rope footpoints (red) located near to 
the active region (middle panel). Later on (right panel), 
reconnection has occurred between the expanding flux 
rope (red), and the neighboring held line (cyan). This 
mechanism migrates the footpoints of the field lines asso- 
ciated with the CME away from the source AR and facili- 
tates a lateral expansion of the CME. The CME footprint 
thus expands to a large extent in the low corona. 
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Figure 8. Panels (la) and (lb) illustrate the scenario 
of the event as proposed by Crooker and Webb [2006]; 
Attnll et al. [2006]; Crooker et al. [2008]. The CME is 
represented by the twisted field line in panel (la) and it 
undergoes ICX with an open field line of the north polar 
coronal hole. Panels (2a) and (2b) show the scenario 
based on the simulation results. The core flux rope first 
reconnects with an overlying closed field line, so closed 
loops (making up the north side lobe), form between the 
polar coronal hole and the north dimming region. At the 
same time, the connectivity of the twisted closed loops 
that constitute the CME is changed, so that a significant 
portion of the CME is effectively displaced and expands 
out into the heliosphere. This extended twisted CME 
field can later reconnect with an open polar field line 
(panel 2c). 



